Unveiling cloudy exoplanets: the influence of cloud model choices on retrieval solutions

In recent years, it has become clear that a substantial fraction of transiting exoplanets have some form of aerosol present in their atmospheres. Transit spectroscopy – mostly of hot Jupiters, but also of some smaller planets – has provided evidence for this, in the form of steep downward slopes from blue to red in the optical part of the spectrum, and muted gas absorption features throughout. Retrieval studies seeking to constrain the composition of exoplanet atmospheres must therefore account for the presence of aerosols. However, clouds and hazes are complex physical phenomena, and the transit spectra that are currently available allow us to constrain only some of their properties. Therefore, representation of aerosols in retrieval models requires that they are described by only a few parameters, and this has been done in a variety of ways within the literature. Here, I investigate a range of parameterisations for exoplanet aerosol and their effects on retrievals from transmission spectra of hot Jupiters HD 189733b and HD 209458b. I find that results qualitatively agree for the cloud/haze itself regardless of the parameterisation used, and indeed using multiple approaches provides a more holistic picture; the retrieved abundance of H2O is also very robust to assumptions about aerosols. I also find strong evidence that aerosol on HD 209458b covers less than half of the terminator region, whilst the picture is less clear for HD 189733b

Exoplanet Transmission Spectroscopy using KMOS

KMOS (K-Band Multi Object Spectrograph) is a novel integral field spectrograph installed in the VLT's ANTU unit. The instrument offers an ability to observe 24 2.8"$\times$2.8" sub-fields positionable within a 7.2' patrol field, each sub-field producing a spectrum with a 14$\times$14-pixel spatial resolution. The main science drivers for KMOS are the study of galaxies, star formation, and molecular clouds, but its ability to simultaneously measure spectra of multiple stars makes KMOS an interesting instrument for exoplanet atmosphere characterization via transmission spectroscopy. We set to test whether transmission spectroscopy is practical with KMOS, and what are the conditions required to achieve the photometric precision needed, based on observations of a partial transit of WASP-19b, and full transits of GJ 1214b and HD 209458b. Our analysis uses the simultaneously observed comparison stars to reduce the effects from instrumental and atmospheric sources, and Gaussian processes to model the residual systematics. We show that KMOS can, in theory, deliver the photometric precision required for transmission spectroscopy. However, this is shown to require a) pre-imaging to ensure accurate centering and b) a very stable night with optimal observing conditions (seeing $\sim$0.8"). Combining these two factors with the need to observe several transits, each with a sufficient out-of-transit baseline (and with the fact that similar or better precision can be reached with telescopes and instruments with smaller pressure,) we conclude that transmission spectroscopy is not the optimal science case to take advantage of the abilities offered by KMOS and VLT.Comment: 11 pages, accepted to MNRA

Constraining the Atmospheric Composition of the Day-Night Terminators of HD 189733b : Atmospheric Retrieval with Aerosols

A number of observations have shown that Rayleigh scattering by aerosols dominates the transmission spectrum of HD 189733b at wavelengths shortward of 1 $\mu$m. In this study, we retrieve a range of aerosol distributions consistent with transmission spectroscopy between 0.3-24 $\mu$m that were recently re-analyzed by Pont et al. (2013). To constrain the particle size and the optical depth of the aerosol layer, we investigate the degeneracies between aerosol composition, temperature, planetary radius, and molecular abundances that prevent unique solutions for transit spectroscopy. Assuming that the aerosol is composed of MgSiO$_3$, we suggest that a vertically uniform aerosol layer over all pressures with a monodisperse particle size smaller than about 0.1 $\mu$m and an optical depth in the range 0.002-0.02 at 1 $\mu$m provides statistically meaningful solutions for the day/night terminator regions of HD 189733b. Generally, we find that a uniform aerosol layer provide adequate fits to the data if the optical depth is less than 0.1 and the particle size is smaller than 0.1 $\mu$m, irrespective of the atmospheric temperature, planetary radius, aerosol composition, and gaseous molecules. Strong constraints on the aerosol properties are provided by spectra at wavelengths shortward of 1 $\mu$m as well as longward of 8 $\mu$m, if the aerosol material has absorption features in this region. We show that these are the optimal wavelengths for quantifying the effects of aerosols, which may guide the design of future space observations. The present investigation indicates that the current data offer sufficient information to constrain some of the aerosol properties of HD189733b, but the chemistry in the terminator regions remains uncertain.Comment: Transferred to ApJ and accepted. 11 pages, 10 figures, 1 tabl

Testing 2D temperature models in Bayesian retrievals of atmospheric properties from hot Jupiter phase curves

Spectroscopic phase curves of transiting hot Jupiters are spectral measurements at multiple orbital phases, giving a set of disc-averaged spectra that probe multiple hemispheres. By fitting model phase curves to observations, we can constrain the atmospheric properties of hot Jupiters such as molecular abundance, aerosol distribution and thermal structure, which offer insights into their dynamics, chemistry, and formation. In this work, we propose a novel 2D temperature scheme consisting of a dayside and a nightside to retrieve information from near-infrared phase curves, and apply the scheme to phase curves of WASP-43b observed by HST/WFC3 and Spitzer/IRAC. In our scheme, temperature is constant on isobars on the nightside and varies with cos$^n$(longitude/$\epsilon$) on isobars on the dayside, where $n$ and $\epsilon$ are free parameters. We fit all orbital phases simultaneously using the radiative transfer package NEMESISPY coupled to a Bayesian inference code. We first validate the performance of our retrieval scheme with synthetic phase curves generated from a GCM, and find our 2D scheme can accurately retrieve the latitudinally-averaged thermal structure and constrain the abundance of H$_2$O and CH$_4$. We then apply our 2D scheme to the observed phase curves of WASP-43b and find: (1) the dayside temperature-pressure profiles do not vary strongly with longitude and are non-inverted; (2) the retrieved nightside temperatures are extremely low, suggesting significant nightside cloud coverage; (3) the H$_2$O volume mixing ratio is constrained to $5.6\times10^{-5}$--$4.0\times10^{-4}$, and we retrieve an upper bound for CH$_4$ at $\sim$10$^{-6}$.Comment: 23 pages, 20 figures, 3 tables, accepted for publication in MNRA

On the potential of the EChO mission to characterise gas giant atmospheres

Space telescopes such as EChO (Exoplanet Characterisation Observatory) and JWST (James Webb Space Telescope) will be important for the future study of extrasolar planet atmospheres. Both of these missions are capable of performing high sensitivity spectroscopic measurements at moderate resolutions in the visible and infrared, which will allow the characterisation of atmospheric properties using primary and secondary transit spectroscopy. We use the NEMESIS radiative transfer and retrieval tool (Irwin et al. 2008, Lee et al. 2012) to explore the potential of the proposed EChO mission to solve the retrieval problem for a range of H2-He planets orbiting different stars. We find that EChO should be capable of retrieving temperature structure to ~200 K precision and detecting H2O, CO2 and CH4 from a single eclipse measurement for a hot Jupiter orbiting a Sun-like star and a hot Neptune orbiting an M star, also providing upper limits on CO and NH3. We provide a table of retrieval precisions for these quantities in each test case. We expect around 30 Jupiter-sized planets to be observable by EChO; hot Neptunes orbiting M dwarfs are rarer, but we anticipate observations of at least one similar planet.Comment: 22 pages, 30 figures, 4 tables. Accepted for publication in MNRA

A uniform analysis of HD209458b Spitzer/IRAC lightcurves with Gaussian process models

We present an analysis of Spitzer/IRAC primary transit and secondary eclipse lightcurves measured for HD209458b, using Gaussian process models to marginalise over the intrapixel sensitivity variations in the 3.6 micron and 4.5 micron channels and the ramp effect in the 5.8 micron and 8.0 micron channels. The main advantage of this approach is that we can account for a broad range of degeneracies between the planet signal and systematics without actually having to specify a deterministic functional form for the latter. Our results do not confirm a previous claim of water absorption in transmission. Instead, our results are more consistent with a featureless transmission spectrum, possibly due to a cloud deck obscuring molecular absorption bands. For the emission data, our values are not consistent with the thermal inversion in the dayside atmosphere that was originally inferred from these data. Instead, we agree with another re-analysis of these same data, which concluded a non-inverted atmosphere provides a better fit. We find that a solar-abundance clear-atmosphere model without a thermal inversion underpredicts the measured emission in the 4.5 micron channel, which may suggest the atmosphere is depleted in carbon monoxide. An acceptable fit to the emission data can be achieved by assuming that the planet radiates as an isothermal blackbody with a temperature of $1484\pm 18$ K.Comment: 18 pages, 5 figures, 6 tables. Accepted by MNRA

Outstanding Challenges of Exoplanet Atmospheric Retrievals

Spectral retrieval has long been a powerful tool for interpreting planetary remote sensing observations. Flexible, parameterised, agnostic models are coupled with inversion algorithms in order to infer atmospheric properties directly from observations, with minimal reliance on physical assumptions. This approach, originally developed for application to Earth satellite data and subsequently observations of other Solar System planets, has been recently and successfully applied to transit, eclipse and phase curve spectra of transiting exoplanets. In this review, we present the current state-of-the-art in terms of our ability to accurately retrieve information about atmospheric chemistry, temperature, clouds and spatial variability; we discuss the limitations of this, both in the available data and modelling strategies used; and we recommend approaches for future improvement

Exoplanet atmospheres with EChO: spectral retrievals using EChOSim

We demonstrate the effectiveness of the Exoplanet Characterisation Observatory mission concept for constraining the atmospheric properties of hot and warm gas giants and super Earths. Synthetic primary and secondary transit spectra for a range of planets are passed through EChOSim (Waldmann & Pascale 2014) to obtain the expected level of noise for different observational scenarios; these are then used as inputs for the NEMESIS atmospheric retrieval code and the retrieved atmospheric properties (temperature structure, composition and cloud properties) compared with the known input values, following the method of Barstow et al. (2013a). To correctly retrieve the temperature structure and composition of the atmosphere to within 2 {\sigma}, we find that we require: a single transit or eclipse of a hot Jupiter orbiting a sun-like (G2) star at 35 pc to constrain the terminator and dayside atmospheres; 20 transits or eclipses of a warm Jupiter orbiting a similar star; 10 transits/eclipses of a hot Neptune orbiting an M dwarf at 6 pc; and 30 transits or eclipses of a GJ1214b-like planet.Comment: 13 pages, 15 figures, 1 table. Accepted by Experimental Astronomy. The final publication will shortly be available at Springer via http://dx.doi.org/10.1007/s10686-014-9397-

A comparison of exoplanet spectroscopic retrieval tools

Over the last several years, spectroscopic observations of transiting exoplanets have begun to uncover information about their atmospheres, including atmospheric composition and indications of the presence of clouds and hazes. Spectral retrieval is the leading technique for interpretation of transmission spectra and is employed by several teams using a variety of forward models and parameter estimation algorithms. However, different model suites have mostly been used in isolation and so it is unknown whether the results from each are comparable. As we approach the launch of the James Webb Space Telescope, we anticipate advances in wavelength coverage, precision, and resolution of transit spectroscopic data, so it is important that the tools that will be used to interpret these information-rich spectra are validated. To this end, we present an intermodel comparison of three retrieval suites: TAUREX, NEMESIS, and CHIMERA. We demonstrate that the forward model spectra are in good agreement (residual deviations on the order of 20-40 ppm), and discuss the results of cross-retrievals among the three tools. Generally, the constraints from the cross-retrievals are consistent with each other and with input values to within 1σ. However, for high precision scenarios with error envelopes of order 30 ppm, subtle differences in the simulated spectra result in discrepancies between the different retrieval suites, and inaccuracies in retrieved values of several σ. This can be considered analogous to substantial systematic/astrophysical noise in a real observation, or errors/omissions in a forward model such as molecular line list incompleteness or missing absorbers